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Abstract 

The distribution of the electron scattering rate on the Fermi surface of a quasi-one-dimensional conductor is calculated 
for the electron-electron umklapp interaction. We find that in certain regions on the Fermi surface the scattering rate 
is anomalously high. The reason for the existence of these "hot spots" is analogous to the appearance of the van Hove 
singularities in the density of states. We employ a generalized r-approximation (where the scattering integral in the Boltzmann 
equation is replaced by the scattering time which depends on the position at the Fermi surface) to study the dependence of 
the electric resistance on the amplitude and the orientation of a magnetic field. We find that the "hot spots" do not produce 
a considerable magnetoresistance or commensurability effects at the "magic angles" . 



In a recent paper [1], Chaikin suggested phenomeno- 
logically that "hot spots" , the regions where the electron 
scattering rate is anomalously high, may be present at the 
Fermi surface (FS) of the (TMTSF) 2 X organic metals. With 
a particular choice of the model parameters, he was able 
to explain the experimental data [2] on the angular depen- 
dence of magnetoresistance (the Lebed's "magic angles" ef- 
fect [3]). However, the microscopic origin of the "hot spots" 
was not specified in Ref. [1]. In an earlier paper [4], it was 
mentioned that the electron-electron umklapp scattering or 
the scattering on a soft phonon mode in the vicinity of a 
charge/spin-density- wave transition may lead to formation 
of the "hot spots" . The latter idea was elaborated in detail 
in Ref. [5]. In the present paper, we calculate straightfor- 
wardly the distribution of the electron scattering rate on the 
Fermi surface (FS) of a quasi-one-dimensional (Q1D) con- 
ductor due to the electron-electron umklapp scattering and 
demonstrate that the "hot spots" indeed exist. 

Let us consider a Q1D metal which consist of parallel 
conducting chains. The axis x is in the direction of the high- 
est conductivity (along the chains), the axes y (also referred 
to as b) and z (or c) are in the directions of the lower and the 
lowest conductivities, and both directions are perpendicular 
to the chains. The Fermi momentum in the ai-direction is 
equal to ±&f (Ji = 1). In the umklapp scattering process, 
the total momentum of the two scattering electrons changes 
by the value 4A?f which is equal to the lattice wave vector in 
the (TMTSF) 2 X compounds. The scattering rate (the in- 
verse of the lifetime r) of the electron with the momentum 
k = (k x , k y , k z ) is given by the following general expression: 
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where gn is the umklapp scattering amplitude, n(e) is the 
Fermi distribution function at the temperature T, and every 
integration is over the Brillouin zone (BZ). 

In a Q1D conductor, the electron dispersion law can be 



written in the linearized with respect to k x form: 

£a(k) = oiVF{k x — «/cf) + 2tb cos(k y + aipi) 
+2t c cos(fe z + aif2). 



(2) 



Here the energy e is counted from the Fermi level, the index 
a = ± labels the electrons whose momenta along the chains 
are close to ±&f, vf is the Fermi velocity, tt and t c are the 
electron hopping integrals between the chains (tt, ^> t c ), the 
transverse momenta k y and k z are dimensionless and vary 
from to 2-7T, and ip\ and <^ 2 are the phases determined by 
the crystal structure of the (TMTSF) 2 X compounds [6]. 

We shall consider only the scattering of the electrons 
located at the FS (e+(k) = 0). In this case, the scattering 
rate becomes a function of only two variables, k y and k z , 
which label positions on the FS; the third momentum k x is 
determined by the condition that (k x , k y ,k z ) belongs to the 
FS. Substituting (2) into (1) and taking several integrals, 
we find: 
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The normalizing coefficient C is introduced in Eqs. (3) and 
(5) in such a way that the function /(£) becomes the Dirac 
function <5(£) when T — > 0. One may conclude that at low 
temperatures the integration over momenta in (4) gives a 
temperature-independent geometrical factor G(k y , k z ), and, 
for all values of (k y ,k z ), 1/r (3) is proportional to T 2 . In 
other words, the scattering rate (3) is a factorized function 
of the temperature and the position at the FS in agreement 
with a standard Fermi liquid theory. However, as we shall 
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Figure 1: Variation of the electron scattering rate 1/r along 
the FS, 2D case. Numbers 1, 2, 3, and 4 correspond to the 
phases tpi = tt/12, tt/4, 0, and n/2; letters a, 6, c, and d 
define the temperature T = tb/240, tb/24, tt/5, and tb- The 
normalizing coefficient 7 is equal to (2n) 4 hv p / 60 gl. 




100 200 300 400 500 
t/T 

b 

Figure 2: Temperature dependences of the electron scatter- 
ing rates for different positions at the FS, 2D case. Squares 
represent the results of the numerical calculations, and the 
solid lines represent the corresponding curve fits. The coef- 
ficient 7 is the same as in Fig. 1. 



see below, the factorization breaks in the two-dimensional 
(2D) case when t c = in Eq. (4). 

Let us consider the 2D case first. In this case, 1/r is 
a function of only one variable k y . Numerically calculated 
according to Eqs. (3)-(6) with t c = and normalized to T 2 , 
the scattering rate l/r(k y ) is shown in Fig. 1 for different 
temperatures T and phases ipi. Since Eq. (4) has the sym- 
metry G(k y ± 7r, k z ± 7r) = G(k y , k z ), the results are shown 
only for a half of the BZ. We observe that, for certain special 
values of k y = k*, 1/t(/c*)T 2 strongly increases as T — > 0. 
In other words, the points k y become increasingly "hotter" 
relative to other points on the FS as temperature decreases. 
Positions of these "hot points" are determined by the phase 
ipi in Eq. (2), namely k^ = ipi. In the exceptional cases 
tpi = and 7r/2, all points on the FS are "hot". 

Temperature dependences of the scattering rate are 
shown in Fig. 2 for different positions at the FS. At the 
positions which are not "hot" , the scattering rate is pro- 
portional to T 2 (curve A); whereas at the "hot points", 
r- 1 ru T 2 \a(t b /T) when tpi ^ ir/4 (curve B) and t _1 ~ 
T 2 (t b /T) 1/3 when Vl = tt/4 (curve C). 

The origin of the "hot points" can be understood in the 
following way. In the 2D case, the geometrical factor G(k y ) 
(4), can be written, in the limit T = 0, as: 

( 2lx ^h m ^ 2) 

G(ky) = J (2n) 2 Z[hF(kv^\k<i\k^)]. (7) 

Integral (7) is divergent for certain values of the parameter 
k y . This divergency is similar to the van Hove singularity 
in the density of states in the 2D case, because Eq. (7) is 
mathematically analogous to the expression for the density 
of states in the 2D case. 

The characteristic function F(k y , <pi] k y x \ ky 2 " 1 ) (6), 
which enters Eq. (7), is shown in Fig. 3. As a function of 
the two variables k^ and k y 2 \ F has six extrema, three of 



them being saddle points. At these points, the first deriva- 
tives of the function F vanish: dF/dky 1 " 1 '^ = 0. Positions 
of the saddle points depend on the values of the parameters 
ky and ip\, and it is always possible to find such a value 
of ky = ky, that the saddle points (SP) appear at the zero 
level of the function F: 

F(*; iVl ;SP) = 0. (8) 

In this case, integral (7) becomes logarithmically divergent 
in the same way as the density of states diverges at the van 
Hove singularity. It turns out that the values of the function 
F at all three saddle points always coincide and are equal 
to: F(k y , (pi; SP) = — 2 sin(2iy5i) sin(fcj, — ipi). Condition (8) 
is satisfied either when ipi = or 7r/2, or when kt = ip\. 
These values are in agreement with the positions of the "hot 
points" found numerically (see Fig. 1). 

When integral (7) is divergent, we should recall that 
originally it was not the Dirac ^-function in this integral, but 
the function / (5) which has a finite width of the order of T. 
The finite width cuts off the logarithmic divergency, so that 
G{ky) ~ \n(tb/T). In the special case ipi = tt/4, all three 
saddle points merge together in the (ky 1 ^, k y 2 ^) space to form 
a "super" -saddle point where not only the first derivatives 
of the function F vanish, but also the second derivatives do: 
d 2 F / dky^ dky^ = 0, i,j = 1,2. In this case, the expansion 
of the function F in the powers of k v in the vicinity of the 
"super" -saddle point starts from the third power, thus the 
temperature cut-off in Eq. (7) produces a stronger depen- 
dence: G(k y ) ~ (tb/T) 1 ^ 3 . These results are in agreement 
with the numerical analysis (see Fig. 2). 

Let us return now to the three-dimensional (3D) case 
when t c ^ in Eq. (4). The numerically calculated scatter- 
ing rate in the 3D case is shown in Fig. 4 as a function of the 
position on the FS. One can notice the warped "hot lines" . 
Their location on the FS is determined by the condition 
that the saddle points of the function tbF(k y , ipi; k y x \ k y 2 ^) + 
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Figure 3: Characteristic function of the scattering rate, 



F(k y ,(fi;ky 1 \ky 2 " 1 ), shown as a function of (k y L> ,k 
k y = cpi = 7r/8. Four BZ cells are shown for clarity. Posi- 
tions of the saddle points are marked by black circles. 
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Figure 4: Variation of the electron scattering rate along 
the FS, 3D case, <pi = ir/4, ip 2 = tt/12, t b jt c = 10, and 
tb/T = 48. The coefficient 7 is the same as in Fig. 1. 
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(see Eq. (4)), defined in the four- 



dimensional (4D) space (A;^ 1 ', k v 2 \ k^\ k z 2 ^), are located at 
the zero level: 



ttsin(2(y3i) sin(fc !/ — 931) + t c sin(2i^2) sin(k z 
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0. (9) 



Saddle points in the 4D space, unlike in the 2D space, are 
integrable, that is they do not produce singularity in the 
function G(k y ,k z ) (4). A non-zero value of t c (t c <C cuts 
off the 2D saddle point divergency, so the ratio between the 
crest and the valley in Fig. 4 is proportional to ln(tf,/t c ) 
(when ip 1 ^ 0,7r/4, and7r/2) or (tb/t c ) 1/s (when 931 = 7r/4). 
The variation of the scattering rate along the "hot lines" 
depends on the phase ip 2 , but never exceeds 15%. At low 
temperatures T < < c , the distribution of the scattering rate 
over the FS stays constant, and r _1 ~ T 2 at all points of 
the BZ. 

Having calculated the scattering rate, let us calculate 
now the magnetoresistance. The linearized stationary ki- 
netic equation in the magnetic field H has the form [7]: 



zHvf dtp 
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(10) 



Here n = no — V> dno /de is the electron distribution function, 
no is the unperturbed Fermi function, the momenta fey and 
k± are parallel and perpendicular to the magnetic field, E 
is the electric field, and v = de/dk is the electron velocity. 
Dependence of the relaxation time r on the variables k\\ and 
k±_ in Eq. (10) reflects the variation of the scattering rate 
along the BZ. Solving Eq. (10) analytically and using the 
result to calculate the electric current, we find the expression 
for the electric resistivity along the chains, p xx ,: 
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(11) 



The orientation of the magnetic field implicitly enters Eq. 
(11) through the definitions of k\\ and k± relative to the BZ. 

The dependence of the resistance on the magnetic field 
H\\c is shown in Fig. 5. It is natural to describe the value of 
the magnetic field using the parameter h which shows how 
many times electrons cross the BZ before being scattered: 



h 



H 2tt 



Ho 



27rc(r- 1 ) 



bevf 



(12) 



where (...) means the averaging over the whole BZ, and 
uj c = bevFH/c is the cyclotron frequency. We observe from 
Fig. 5 that the resistance p xx (H) saturates at h > 1, when 
electrons cross the whole BZ at least once during their life- 
time. The following relation between the limiting values of 
the resistance can be derived from Eq. (11): 



p xx (H = 00) / Pxx (H = 0) = <r> / ((r- 1 )^)*,. 



(13) 



where (...)jb|| ± is the average taken in the direction parallel 
(||) or perpendicular (_L) to the magnetic field over the BZ. 
In the 2D case, when H\\c, Eq. (13) takes a simpler form: 



p xx (H = 00) /p xx (H = 0) = (t) (1/t). 



(14) 



Due to the small variation of the scattering rate along the 
k z axis (see Fig. 4), the saturation value of the magnetore- 
sistance in Fig. 5 is in agreement with the 2D formula (14). 
The calculated magnetoresistance is too small (13%) com- 
pared to the experiment [2] which shows a much higher value 
and no signs of saturation in a high magnetic field. 

Since in the 3D case the scattering rate has the uniform 
temperature dependence T 2 in the low temperature region 
(T < t c ), the ratio p xx (H = 00) /p xx (H = 0) (14) docs not 
depend on temperature in this region. However, in the 2D 
case, the different temperature dependence of r _1 (&*) at 
the "hot points" can produce an increase of the magnetore- 
sistance as T — > 0. In Fig. 6, we show the maximum value 
of magnetoresistance as a function of the inverse tempera- 
ture for the different phases ipi in the 2D case. We observe 
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Figure 5: Magnetoresistance as a function of the magnetic 
field H\\c, 3D case, Ho = 27rc(r _1 )/6e«F- 



Figure 6: Temperature dependence of the saturation value 
of the magnetoresistance, 2D case. 



that even in the "super" -saddle case (cpi = it/ 4), the tem- 
perature dependence of the magnetoresistance is weak and 
practically negligible. 

In Fig. 7, the dependence of the magnetoresistance on 
the inclination angle of the magnetic field, 8, is shown in the 
situation when the magnetic field is rotated perpendicular 
to the chains. In agreement with the experiment [2], Fig. 
7 shows no magnetoresistance when H\\c, because in this 
orientation the Lorentz force does not sweep electrons into 
the "hot lines". However, in contradiction with the experi- 
ment, Fig. 7 shows no dips or any features at magic angles 
[3] when tan 8 is equal to a rational number. The curves 
in Fig. 7 turn out to be smooth because the scattering rate 
varies too weakly in the k z direction in Fig. 4. 

We conclude that a strong variation of the electron scat- 
tering rate along the Fermi surface should exist in Q1D con- 
ductors in the case of electron-electron scattering. The rea- 
son for the variation is very general and analogous to the 
reason for the existence of the van Hove singularities in the 
density of states. However, we found that the variation is 
not strong enough to produce the high magnetoresistance or 
the angular commensurability effects observed experimen- 
tally. In our consideration we neglected the influence of the 
magnetic field on the scattering rate. That is correct in the 
classical region T ^> uj c . However, at low temperatures the 
influence, which is a quantum effect, must be taken into 
account as it was done in Ref. [3]. We conclude that the 
classical effect of the magnetic field, that is, a specific or- 
der of averaging of the scattering rate due to the Lorentz 
orbital motion of the electrons, is not enough to explain the 
peculiarities of the transport properties of the (TMTSF^X 
conductors in a high magnetic field. 
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Figure 7: Angular dependence of the magnetoresistance for 
different values of the magnetic field, 3D case. The defini- 
tion of h is given in Eq. (12). 
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